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ABSTRACT 

Single atoms can be considered as the most basic 
objects for electron microscopy to test the micro- 
scope performance and basic concepts for modeling 
image contrast. In this work high-resolution trans- 
mission electron microscopy was applied to image 
single platinum, molybdenum and titanium atoms in 
an aberration-corrected transmission electron micro- 
scope. The atoms are deposited on a self-assembled 
monolayer substrate which induces only negligible 
contrast. Single-atom contrast simulations were 
performed on the basis of Weickenmeier-Kohl and 
Doyle- Turner form factors. Experimental and simu- 
lated image intensities are in quantitative agreement 
on an absolute intensity scale which is provided by 
the vacuum image intensity. This demonstrates that 
direct testing of basic properties like form factors be- 
comes feasible. 



application fields like catalysis and crystal growth. 
Apart from application relevance, single atoms can 
be considered as the most basic objects in electron 
microscopy. Imaging single atoms by electron mi- 
croscopy has been previously demonstrated by scan- 
ning transmission electron microscopy (STEM) for 
some t ime already dBatson et al.l . l2002t ICrewe et all 



1970; iKrivanek et all l2010h . Single atoms 



are 

visualized in the annular dark-field STEM mode 
due to the small electron probe formed by the 
condenser-lens system and the atomic-number de- 
pendent scattering-angle distribution. 

Alternatively, high-resolution transmission elec- 
tron microscopy (HRTE M) can be applied for s ingle- 
atom imaging (lliiimal . Il977t iMever et all . 120081: 
lOhnishi et all Il998l) HRTEM is based on phase 
contrast where the interaction between electron 
wave and atom only induces a phase shift with 
negligible amplitude variation in the image elec- 
tron wave. This leads to weak single-atom contrast 
which depends strongly on the defocus and other 
lens aberrations. Comparison between experimen- 
tal and simulated images taken at different defocus 
values is therefore indispensable for the interpreta- 
tion of HR TEM images. This was shown by sev- 
eral groups ( Hi final . Il977t iKoizumi et all l200ll ) who 
compare characteristic features of experimental and 
simulated single-atom images or relative intensities 
between atoms in an image. 

However, the comparison must be performed on 
an absolute intensity scale to ultimately analyze mi- 
croscope performance, basic concepts for modeling 
image formation or even form factors. Absolute in- 
tensity scaling can only be achieved by image nor- 
malization with the vacuum intensity rather than in- 
tensity normalization with respect to another sam- 
ple feature. In this sense, we present for the first 
time quantitative agreement between a series of ex- 
perimental and calculated single-atom images on an 
absolute intensity scale. The agreement is found 
for single exposure images, and no averaging over 
several images is needed to enhance the single- atom 
contrast. Our work is facilitated by an correction- 
lens system for spherical and other aberrations (for 
short Cs-corrector) which improves the resolution to 
better than 0.1 nm and allows to measure precisely 
residual lens aberrations. 



INTRODUCTION 

Imaging single atoms has always been a challeng- 
ing goal in electron microscopy which is interesting 
with respect to studying single-atom dynamics in 
context with cluster nucleation and coarsening in 



MATERIALS AND METHODS 

The samples were prepared by electron-beam 
evaporation of approximately 0.2 monolayers Pt (1 
monolayer: 1.3 • 10 15 atoms /cm 2 ), Mo or Ti, on 
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TABLE I. Measured residual aberrations for the Pt- 
specimen images (Z: defocus, C3, Cs: spherical aber- 
ration, Al: twofold astigmatism, A2: threefold astigma- 
tism, A3: fourfold astigmatism, A4: fivefold astigma- 
tism, B2: axial coma, S3: star aberration). 
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FIG. 1. Modulation transfer function of the Gatan Ultra- 
scan 1000 CCD-camera determined by the edge method. 



a self-assembled monolayer (SAM) film denoted as 
nanosheet in the following. It consists of 1,1'- 
biphenyl-4-thiol molecules (C12H10S) with a thick- 
ness of about 1.0 nm which are crosslink ed by elec- 
tron irradiation (jTurchanin et all l2009h . HRTEM 
was carried out with an aberration-corrected FEI 
Titan 3 80-300 microscope operated at 300 keV. A 
defocus series was acquired which consists of 36 sin- 
gle images with a defocus step of 1 nm between 
19 nm overfocus and —16 nm underfocus. The ex- 
posure time was 1 s. The Cs-corrector was used to 
minimize and determine the residual objective-lens 
aberrations on the basis of t he procedure described 
by U hlemann and Haider ( Uhlemann and Haideii 
Il998l ). Tab U shows the aberration values which were 
used for the image simulations in the case of the Pt- 
atom images. 

Image simulations are based on an object exit 
wave fun ction t/>q calculated with th e multislice for- 
malism (ICowlev and Moodid. Il957l) of the ST EM- 
sim program (jRosenauer and Schowalterl [2006). A 
region around a single Pt-atom was defined large 
enough to avoid artifacts that arise from the pe- 
riodic boundary conditions implied in the mul- 
tislice algorithm. The atomic potential is cal- 



culated by inverse Fourier transformation from 
Weickenmeier-K ohl dWeickenmeier and Kohl I1991D 
or Doyle- Turner (jDovle and Turned . 11968) form fac- 
tors. The former explicitly deals with absorption 
arising from inelastic scattering. The aberrations 
of the imaging-lens system of the microscope mod- 
ify the phase of ipo as a function of the spatial 
frequency g which is described by the wave trans- 
fer function T(g). The image wave function ipi 
is given in Fourier representation by ipi = ipo ' 
T(g) = ipo ■ exp(— ix(g))with the spatial frequency- 
dependent phase shift \(g) elaborated up to t h e fifth 
order by Uhlemann (jUhlemann and Haideii Il998l) 
and Kirkland (?). The partial spatial coherence of 
the incident electron wave is taken into account by 
multiplying T(q) with a damping envelope function 
( Wade and Frankill977D . which is negligible for the 
small Cs and Z values used in this work. Partial 
temporal coherence is characterized by the param- 
eter defocus spread A which is ss 4.5 nm for our 
electron microscope. The effect of partial temporal 
coherence was taken into account by simulating im- 
ages for nineteen different defocus values around the 
chosen Z value and weighting the different contribu- 
tions to the image according to a Gaussian distribu- 
tion with a FWHM given by A. The FEI Titan 3 80- 
300 is equipped with a monochromator which, how- 
ever, was not excited in the present work. Slight il- 
lumination inhomogeneities may occur due to an ad- 
ditional solid aperture implemented in the gun sec- 
tion of the microscope. Image detection by the CCD 
camera dampens pixel-size contrast like a low-pass 
filter. This behavior is described by the modulation 
transfer function (MTF) whic h was measured by 
the ed ge method described by IWeickenmeier et al.l 
(1995) and taken into account by convoluting the 
image with the MTF. The MTF for the Gatan Ul- 
trascan 1000 CCD camera is shown in FigJT] It was 
determined from images taken under the same ex- 
perimental conditions than the images presented in 
this work. Nyquist frequency corresponds in the case 
of Pt-atom images to « 20 nm -1 . 



RESULTS 

Fig. [2] shows typical HRTEM images of small pla- 
nar Pt-clusters and isolated Pt-atoms on a nanosheet 
at different Z-values as indicated on the images. The 
spatial sampling rate was 0.0245 nm/pixel. Fig. [2^) 
shows bright atom contrast under overfocus condi- 
tions. The intensity of the Pt-atom contrast vanishes 
almost completely at Z = 5 nm (Fig. HJd)), which 
is the image with the lowest contrast that could be 
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quantified. At a defocus value of Z = 2 — 3 nm 
no contrast is observed. The minimum is shifted 
into the overfocus region due to the finite positive 
value of Cs- Contrast is reversed under underfo- 
cus conditions (Fig. [2}:,d)) where dark atom con- 
trast is observed. The atoms are highly mobile 
on the nanosheet which can be inferred from the 
change in the atom configuration of the clusters. 
Fig. W?) shows a magnified region where isolated 
atoms can be discriminated from atoms in small 
clusters. Furthermore a comparably large defocus 
of Z = — 19 nm allows distinguishing the contrast 
associated with the nanosheet from the vacuum as 
shown in Fig. [^Jj). The top right corner shows a 
small hole in the nanosheet marked by an arrow. 
Fig. [5f) presents intensity line profiles measured 
at Z = — 8 nm along a line across a vacuum and 
nanosheet region with mean intensities correspond- 
ing to 2994 and 3038 counts. Almost the same 
standard deviation of the intensity (90 counts for 
vacuum and 110 counts for the substrate) for these 
two regions indicates the uniform small thickness of 
the nanosheet and its negligible influence on image 
contrast at small defocus values. The analysis of 
the intensity profiles yields a noise level of about 
3 %. Therefore, the nanosheet contrast can read- 
ily be neglected in the simulations for the defocus 
range observed. Normalization of the experimental 
images with the nanosheet background intensity cor- 
responds in a good approximation to normalization 
with the vacuum intensity. Nonetheless, a region 
in direct vicinity of the analyzed atom should be 
taken for normalization. Due to the monochroma- 
tor of the microscope, the illumination is not uni- 
form over the field of view which can be seen from 
the gradient in the images in Fig. [2^,)-d). For several 
different defocus values intensity profiles of isolated 
atoms were obtained from experimental images by 
averaging profiles along the horizontal and vertical 
image axis. Fig. [3] shows representative experimen- 
tal intensity profiles at different defocus values. In- 
tensity profiles of simulated single-atom images were 
extracted in the same way and are also contained in 
Fig. [3J Experimental and simulated profiles agree 
well with respect to the noise level of 3 %. For quan- 
titative comparison of experimental and simulated 
images, the maximum intensity normalized with re- 
spect to the background intensity and the full width 
at half maximum (FWHM) of the peak were chosen. 
Another quantitative feature is the radius of the first 
side minimum/maximum (overfocus/underfocus im- 
ages) of the fringe-like contrast around the central 
peak. 

Simulated data for the peak intensity values, 



FWHM and distance (radius) from the center 
to the first side minimum/maximum (underfo- 
cus/overfocus image) are listed in Table [TTI together 
with the experimental data for twelve different Z val- 
ues. Within the signal-to-noise ratio full agreement 
between simulated and experimental peak intensities 
is obtained for Weickenmeier-Kohl form factors. For 
Doyle- Turner form factors, a different defocus de- 
pendence is observed as shown in the forth column 
of Table El 

For Weickenmeier-Kohl form factors good agree- 
ment is also found for the radius of the first side 
minimum/maximum as a result of the precise de- 
termination of objective lens aberrations. Quantita- 
tive agreement between simulation and experiment 
is also demonstrated in Fig. [3] It shows six represen- 
tative intensity profiles for three different overfocus 
values in Fig. [3^,) and three different underfocus val- 
ues in Fig. [5]d). Peak intensity and FWHM are in 
good agreement. Overfocus profiles clearly show the 
fringe-like contrast around the central maximum. 
We emphasize that quantitative agreement between 
experimental and simulated single-atom contrast is 
achieved on an absolute scale only if the MTF of 
the CCD camera is taken into account. The impor- 
tance of the MTF was also ou tlined in cont ext with 
the "Stobbs-factor" problem (jThustl [2009]) . where 
quantitative agreement for crystalline specimens is 
found when accouting for the MTF. 

The method for quantitative single-atom imag- 
ing shown in this report is generally applicable to 
all elements which can be deposited in very small 
amounts onto a nanosheet by electron-beam evapo- 
ration. Apart from platinum (atomic number: 78) 
we also analyzed molybdenum (atomic number: 42) 
and titanium (atomic number: 22). Mo is a signif- 
icantly weaker scatterer than Pt resulting in a re- 
duced atom contrast which is more subject to noise. 
Nonetheless a quantitative comparison between ex- 
perimental and simulated single-atom intensity pro- 
files is possible. A defocus series of the Mo-sample 
was taken in the same way as for Pt apart from 
magnification that was increased to achieve a sam- 
pling rate of 0.0138 nm/ pixel. Due to the strong 
effect of the MTF, the contrast of Mo-atoms cannot 
be directly compared to the contrast of Pt-atoms, 
and should only be compared with image simulations 
performed under the same conditions (i.e. aberra- 
tions, magnification). From the experimental defo- 
cus series intensity profiles for isolated single Mo- 
atoms were taken in analogy to Figj3] for a defocus 
range from Z = —10 nm to Z = +10 nm. The pro- 
files were then stacked in y-direction into an image 
which is shown in Figure |4j The left part of the im- 
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TABLE II. Comparison of contrast characteristics between experimental and simulated images of single Pt-atoms 
at 12 different defocus values (1st column). Columns 2-4: peak contrast of single Pt-atoms with respect to the 
background intensity in percent for experimental und simulated data for Weickenmeier-Kohl (WK) and Doyle- Turner 
(DT)form factors. Column 5-8: FWHM of central peak and radius of the first side minimum (overfocus images) and 
first side maximum (underfocus images) ring in experiment and simulation. The error for distance measurements is 
approximately 0.013 nm corresponding to half a pixel. 
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age shows the experimental profiles. The right side 
shows profiles from simulated defocus-series images. 
The missing defocus values in the experimental data 
are due to the low contrast of Mo-atoms in these 
images. Apart from defocus value Z = 6 nm atom 
contrast is in good qualitative and even quantitative 
agreement. In the underfocus regime, the bright halo 
around the dark atom contrast is reproduced as ob- 
served in the simulation. This feature corresponds 
to the first side maximum already discussed earlier. 

In the case of Ti the comparison becomes more dif- 
ficult mainly due to the nanosheet contrast, which 
hampers the possibility to find isolated atoms and 
to normalize the contrast with respect to vacuum in- 
tensity. Figure [5] shows Ti-atoms on the nanosheet 
at a rather high defocus of Z = —10 nm. Due to 
the comparably weak scattering power of the Ti- 
atoms, the Ti-contrast does not significantly exceed 
the contrast of the nanosheet anymore. This pre- 
vents normalization of the atom contrast with re- 
spect to the background intensity. For smaller de- 
focus values single-atom contrast is even further di- 
minished but the noise remains the same. It should 
also be noted that the binding to other atoms can 
slightly change ima ge contrast for lighter elements 
( Meyer et all 1201 lh . Based on the quantification of 
Mo-atom contrast and the comparison to Ti-atoms 
we conclude that the lightest elements, which al- 
low quantification of single-atom contrast, should 
be elements in the range of atomic number 30. 
This is mainly due to the increasing contrast of 



the nanosheet with respect to the single- atom con- 
trast, but also due to the reduction of the signal- 
to-noise ratio which can be achieved in a single- 
exposure image of light single atoms. Furthermore 
we note, that each biphenyl molecule contains a S- 
atom (atomic number: 16) which cannot be distin- 
guished from Ti-atoms (atomic number: 22). How- 
ever the nanosheet may contain only a small concen- 
tration of S -atoms due to the nano sheet preparation 
procedure ( Turchanin et all 120091) . 



DISCUSSION 

In the previous section, quantitative agreement 
between experimental and simulated single-atom 
contrast was shown. Negligible contrast of the sub- 
strate is essential for normalization of single-atom 
contrast with respect to vacuum intensity. This 
is a prerequisite to quantitatively compare experi- 
mental and simulated single-atom intensities on an 
absolute scale and to test form factors which are 
used in all the image simulation software packages. 
Weickenmeier-Kohl (WK) and Doyle- Turner (DT) 
form factors are established for high-energy electron 
diffraction and are therefore tested in this work for 
wave-function simulations. Elastic scattering is de- 
scribed by a real form factor f e i while an imaginary 
(absorptive) part f a bs was introduced by Weicken- 
meier and Kohl (1991) to take into account ther- 
mal diffuse scattering in crystals. A better agree- 
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ment between simulated and experimental images 
of single Pt-atoms is obtained for the WK form fac- 
tor which is surprising on first sight because ther- 
mal diffuse scattering should be irrelevant for sin- 
gle atoms. This peculiarity can be understood by 
considering in detail the calculation of / e ; for the 
WK and DT form factors. Although the physics be- 
hind f e i is identical, the fit functions for / e ; differ 
for the WK and DT form factors. The elastic part 
of the WK form factors is calculated on the basis of 
fei(s) = s~ 2 £\ =1 Ai[l - exp(-BiS 2 )] as a function 
of s = g/Air with parameters Aj and Bi given by 
Weickenmeier and Kohl (1991). Doyle and Turner 
(1968) suggested f el (s) = £) i=1 Aiexp(-BiS 2 ) with 
different parameters A and Bi which are listed, 
e.g., in the Inte rnational Tables for Crystallography 
(ICowlevl . fl99l . The different fit functions and pa- 
rameters Ai , Bi lead to higher scattering amplitudes 
at small spatial frequencies for f e i,WK which is the 
main reason for the good agreement of the WK form 
factors with the experimental data. The influence of 
the imaginary part of the WK form factor is negli- 
gible for single-atom scattering. In both cases, the 
form factors are dampened by a Debye- Waller expo- 
nent exp(—Mg 2 ) which contains the Debye Waller 
factor M. However, t he variation of M b e tween , e.g., 
A 2 and 0.046 A 2 (|Sears and Shellevt Il99lh only 
leads to a negligible reduction of the amplitude of the 
scattered electron wave for a single Pt-atom. Since 
the argument of the exponential function is propor- 
tional to — g 2 the influence of M is negligible for the 
small g-values which are important for single-atom 
imaging. 

While form factors are the basis for calculating 
electron-wave functions of an atom, image contrast 
depends also on instrumental parameters. These 
are magnification, objective lens aberrations, par- 
tial incoherencies as well as the modulation transfer 
function of the imaging detector. All these param- 
eters must be known as precisely as possible during 
the experiment for quantitative studies. These are 
then used in the subsequent image simulation pro- 
cedure. The Cs-corrector reduces spherical aberra- 
tions to a point where the main influence on im- 
age contrast is defocus. Variation of defocus allows 
analysis of the defocus dependence of image con- 
trast. If image aquisition is performed with a CCD- 
camera, the MTF has a strong effect on image con- 
trast, and quantitative agreement is only obtained if 
the MTF is determined for the camera and camera 
mode which is used for the experiment. Finally, par- 
tial coherence also has an effect on image contrast. 
Partial spatial coherence can be neglected, because 



primary objective lens aberrations are small due to 
the Cg-corrector. Partial temporal coherence was 
calculated from instabilities of high tension and lens 
current, which turned out to be precise enough for 
an agreement of image contrast. 



SUMMARY 

To summarize, through-focal series of HRTEM 
images of single Pt-, Mo- and Ti-atoms were ac- 
quired with an aberration-corrected transmission 
electron microscope. The samples were prepared by 
electron-beam evaporation on self-assembled mono- 
layer substrates wich provide negligible background 
contrast as a prerequisite for quantitative analysis of 
single-atom contrast on an absolute intensity scale. 
The imaging process is modeled by Fourier optics 
taking into account all known instrumental effects 
including wave aberrations, partial coherence and 
the modulation transfer function of the CCD cam- 
era. Doyle- Turner and Weickenmeier-Kohl form fac- 
tors are compared for wave-function simulations. We 
find quantitative agreement between simulated and 
experimental single-exposure image intensity of sin- 
gle Pt- and Mo- atoms based on Weickenmeier-Kohl 
form factors which provide a more adequate mod- 
eling of the real part of the form factor compared 
to the Doyle- Turner form factors. The influence of 
Debye- Waller dampening and the absorptive part of 
the Weickenmeier-Kohl form factors is negligible for 
single-atom images. We emphasize that the agree- 
ment is found on an absolute intensity scale based on 
raw image data without manipulation of the images 
like filtering scaling or brackground substraction. 
Single-atom contrast of Ti-atoms with an atomic 
number of 22 is not sufficiently strong with respect 
to the nanosheet contrast. We estimate that the pre- 
sented method allows single-atom contrast quantifi- 
cation down to an atomic number of approximately 
30, the limiting factor mainly being signal-to-noise 
ratio. 
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FIG. 2. Experimental HRTEM images and substrate contrast. a)-d): HRTEM images of single Pt-atoms and small 
clusters of Pt-atoms on a nanosheet acquired with a FEI Titan 3 80-300 microscope at 300 kV and different defocus 
values Z as indicated, e) Magnified and strongly defocused (Z = — 19 nm) image. Single Pt-atoms and substrate 
are clearly visible, f) Intensity profiles taken from a HRTEM image d) taken at Z = —8 nm in a vacuum (gray) and 
nanosheet (black) region. Intensity (absolute counts) and noise levels are comparable. 
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FIG. 3. Comparison of simulated and experimental 
intensity profiles. Experimental (gray) and simulated 
(black) intensity profiles across single Pt-atoms at a) 
overfocus values of 7 nm, 6 nm and 5 nm and b) un- 
derfocus values of -1 nm, -2 nm and -3 nm. 



0.138 -0.138 
position [nm] 



FIG. 4. Comparison of experimental and simulated in- 
tensity profiles of single Mo-atoms as a function of defo- 
cus. Gray-scale values of the images correspond to abso- 
lute image intensities. The left side shows horizontally 
intensity profiles stacked vertically for different dcfocus 
values between Z = 10 nm and Z = — 10 nm. The miss- 
ing part at defocus values around Z — nm is due to 
the low contrast of Mo-atoms. The right side of the im- 
age shows the corresponding simulated intensity profiles. 
Good agreement is achieved except for Z = — 6 nm. 
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FIG. 5. Experimental HRTEM image of deposited Ti- 
atoms on a nanosheet, imaged with a defocus value of 
Z — — 10 nm. Isolated Ti-atoms are clearly visible due to 
dark atom contrast. The contrast of nanosheet is strong 
with respect to vacuum intensity and can therefore not 
be neglected. 
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